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Abstract 

If the electroweak symmetry breaking is originated from a strongly coupled sector, 
as for instance in composite Higgs models, the Higgs boson couplings can deviate 
from their Standard Model values. In such cases, at sufficiently high energies 
there could occur an onset of multiple Higgs boson and longitudinally polarised 
electroweak gauge boson (Vl) production. We study the sensitivity to anomalous 
Higgs couplings in inelastic processes with 3 and 4 particles (either Higgs bosons 
or Vls) in the final state. We show that, due to the more severe cancellations 
in the corresponding amplitudes as compared to the usual 2 — > 2 processes, 
large enhancements with respect to the Standard Model can arise even for small 
modifications of the Higgs couplings. In particular, we find that triple Higgs 
production provides the best multiparticle channel to look for these deviations. 
We briefly explore the consequences of multiparticle production at the LHC. 
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1 Introduction 



The search for the mechanism responsible for electroweak symmetry breaking (EWSB) has 
been a long-lasting endeavour that may finally be resolved with the new data from the 
LHC However, even with the recent evidence for a light scalar state, the jury is still out on 
whether EWSB is caused by new strong interactions or not. In order to have a definitive 
answer further experimental probes are necessary. 

In particular, a hallmark of strong interactions in the EWSB sector is multiple particle 
production [l] [j] In a strongly coupled EWSB sector one would expect copious production of 
longitudinal gauge bosons granted enough energy is available to produce them. This would 
be similar to the production of events with a large pion multiplicity in QCD at high energies. 
In fact, multi-W production was studied in a simplified scaled-up version of QCD almost 20 
years ago (3]. In this Letter we study the inelastic production of longitudinally polarised W 
and Z bosons (denoted collectively by V£) and Higgs bosons in the context of an effective 
Lagrangian. We estimate the energy scale at which these processes become relevant, which 
signals the onset of new physics, as recently discussed in |4). In particular, we will be inter- 
ested in the sensitivity to non-SM Higgs couplings in the growth of the cross section for these 
processes. Our results have as a particular case the study of unitarity violation in multi-Vi 
production in the Higgsless model [5). We show that models with partial unitarization, such 
as the composite Higgs model, can lead to a large enhancement of multiparticle cross section 
due to the absence of cancellation mechanisms in the corresponding scattering amplitudes. 
This effect becomes more acute as the final state multiplicity increases, provided that enough 
energy is available. 

2 Multiparticle cross sections and unitarity 

The perturbative unitarity bound in the inelastic 2 — » n process assuming s-wave dominance 
for a given center-of-mass energy yfs is |5}m: 



This unitarity bound sets stringent constraints on the scattering amplitudes. Since the 
relativistic n— body phase space is proportional to s n ~ 2 the unitarity bound requires that 
the amplitude grows with energy no faster than 



One should note, however, that multiple particle production with large cross section could also be 
obtained in weakly coupled theories at tree level simply due to the large number of diagrams, but this is 
expected only for very large multiplicities, of order 0{\/oiew) PI- 





s 




(2) 
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Suppose, for instance, that the EWSB sector is described by a simple nonlinear sigma model 
(NLaM), neglecting transverse gauge bosons for the moment, and assuming possible reso- 
nance states to be very heavy: 

C NLaM = ^Tr [drU&rf] (3) 
where v = 246 GeV is the usual scale of electroweak symmetry breaking and 

U = e—. (4) 



The isospin triplet "pion" fields 7r* {% = 1, 2, 3) will be identified with the longitudinally po- 
larised gauge bosons through the equivalence theorem [l] . By power-counting, the scattering 
amplitude in this model grows with energy as 

A NLaM (2 -»■ n) ~ 4" ( 5 ) 
v 

and hence naively 

*( 2 ^)~l(^) 2sn ~ 2 - ( fi ) 
Therefore, the growth of the cross section towards the unitarity bound in this model is faster 

for larger number of particles due to the kinematical factors in the phase space, assuming of 
course that enough energy is available. 

Conversely, there must be stronger cancellations in the scattering amplitudes due to new 
physics as the number of final state particles is increased. For instance, unitarity requires 
that A(2 — > 2) ~ constant and A(2 — > 4) ~ 1/s, whereas they both grow as ~ s in the 
NLaM. Therefore, in the absence of a perfect cancellation, it is plausible that the growth of 
the cross section may have a large impact in multi-V^ production. It is the purpose of this 
work to examine this impact. 



Given the fully relativistic n— body phase space given by [7]: 



*•<«> = / n j^m-s^mvi - ± Pl) = g9|^g s- <7) 



one can easily estimate the energy scale A n at which perturbative unitarity is violated in 
2 — > n processes in the NLaM: 



A, 



2(n - l)\(n 



In 

V. 



_(27r) 3 - 3 "(7r/2)"- 

For example, unitarity is violated in 2 — > 4 processes at an energy which is almost 2.5 higher 
than that for the usual 2—7-2 processes. This estimate is in reasonable agreement with (5]. 
One should notice that we are not including in this rough estimate the growth due to the 
combinatorial factors and a proper phase space integration. These will be included below in 
a fully automated calculation. 
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3 Anomalous Higgs couplings and partial unitarization 



In order to recover unitarity there must an UV completion of the model describing the 
interaction of the lightest degrees of freedom. The simplest possibility is the addition of a 
scalar field, which is identified with the Higgs scalar. However, it is possible that the Higgs 
scalar is a composite particle with couplings that may differ from the SM ones. In this case, 
the theory is not UV-complete and unitarity is only partially restored. For such a theory 
one can use an effective Lagrangian (SMEFF) to parameterise its couplings to longitudinally 
polarised gauge bosons and self-couplings (couplings to fermions are not relevant to the 
results presented here) (8): 



eff 



+ 



1 + 2a- 



h 

i 

v 
1 



h 2 , h 3 
+ h— + 



•v 



Tr [d^Ud^tf] 



2 l2 

m h h 



A 



d 3 Xvh 3 - d 4 -h 4 + 



(9) 



2 y p ' 2 

This parameterisation is common to a large class of models, such as composite Higgs models, 
and has been used to study anomalous Higgs couplings in VlVl —> VlVl, hh processes at the 
LHC (8 ,11, 19 . Unitarity is recovered for the SM values a = b = ds = d 4 = l and 63 = 0. 
For different values of these parameters the usual cancellation provided by the scalar field is 
incomplete. As an example, the Minimal Composite Higgs Model (MCHM4) predicts that 
the couplings of the "pions" with the Higgs boson follows from an expansion around the 
vacuum h(x) = of the effective Lagrangian PL?] 



P 



sin 



e 



h(x) 

~T 



Tr [d^Ud^tf] 



(10) 



with the identification v = f sin 9, which comes from the mass term for the gauge fields. This 
Lagrangian has a discrete symmetry under the parity transformation h — > —h and n — > —71, 



although this is not obvious in this representation 12 . Therefore, in the MCHM4 
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v/T^C; 6 = l-2£; b 3 



In order to study the 2 — > 4 scattering, one must expand each field U to order 0{n & ): 



15w 2 



-7T ■ 7T 



(12) 



The number of diagrams increases considerably with the number of final state particles, 
making it impractical to perform an analytical computation. Hence we have implemented 
the Lagrangian given in Eq.g in FormCalc (13) and MadGraph (14) using FeynRules (15 
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(with UFO output 16 for the higher dimensional operators). We have also implemented 
this model in CalcHEP 17 using LanHEP package 18 with the help of auxiliary fields. 
This Lagrangian is equivalent to the usual linear sigma model by a field redefinition. 



In the familiar case of 2 — > 2 amplitudes the only kinematical dependence is on the 
Mandelstam variables s and t. For instance, denoting the goldstone bosons by their electric 
charge, the 00 — > H — amplitude, arising from only 4 diagrams, is given by: 



Hence one can easily see that there is a violation of unitarity even with the presence of the 
Higgs boson if its coupling is not SM-like, i.e., a ^ 1. However, in the SM one obtains a 
constant amplitude at high energies, as expected. 

The 2 — > 4 amplitudes are much more complicated, containing of the order of 100 di- 
agrams and depending on several combinations of the scalar products of the different 4- 
momenta involved. However, some of their properties can be demonstrated with the follow- 
ing simple example for a given point in phase space, where all the particles are in the same 
plane (we will keep the "pions" massless at the amplitude level since their masses are not 
relevant for issues of unitarity), in which case we obtain 



It grows with s, as expected. However, in the SM (a = b = d% = 1) one obtains in the limit 
s > ml: 



and we explicitly see the strong cancellation where the first two terms in the amplitude 
vanish and the behaviour change from s/v A to m^/(st> 4 ), as anticipated. The triple Higgs 
anomalous coupling parameterised by d% does not enter in the dominant contribution. In the 
following we will take d 3 = 1. There is no contribution from the couplings d^ and 63 for the 
above processes. This result for the amplitude depends on the phase space configuration. 
The polynomials will be different but the features described above also happens at other 
phase space configurations and other channels. 

For 2 — y 3 processes a similar analysis can be performed. For instance, again for a given 




(13) 



AVoo+- oc £ [72s (13a 4 - a 2 (7b + 5) - 1) + 
3ml (1580a 4 - 378a 3 d 3 - 3a 2 (2456 + 131) - 74) + 
^ (9774a 4 - 3087a 3 d 3 - a 2 (44946 + 1289) + 52) + 



(14) 



- A/ f 0;00+- (X 




(15) 
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configuration in phase space we find 

M 00 ;hhh oc ^3 [s (-4a 3 + Aab - 36 3 )) - 
m\ (-8a 3 + 8a6 + 36 3 ) + 
^(a 3 + a6- 66 3 - 3a 2 d 3 ) + ---l (16) 

and 

Moo.+- h oc -qf-3 [a (-1 + 2a 2 - 6) + 



192v d 

^ (-164 + 386a 2 - 2136 - 9ad 3j 



^ (-262 + 291a 2 - 936 + 81ad 3 ) + • • • (17) 

We again find that for the SM the first two terms in these amplitudes vanish, as it should. 
Notice also that the Aioo-,hhh amplitude is sensitive to 6 3 , being the lowest multiplicity process 
in which this happens. In addition, as can be anticipated from the parity of the MCHM4 
class of theories, under which n —■ —n and h — > —h [l2j, the polynomial with largest growth 
in the 2—7-3 processes also vanish for the values of a, b and o 3 that obey the MCHM4 



relations as in Eq.(ll) 



In summary, for a ^ 1 and b ^ 1 the squared amplitude grows as s 2 instead of decreasing 
as 1/s 2 . Therefore there is a large sensitivity of 2 — > 3 and 2—7-4 processes to non-SM Higgs 
couplings compared to 2 — > 2 processes, whereas the SM amplitude goes to a constant for 
large s. For the 2—7-3 processes, there is a also a suppression if the values predicted by the 
MCHM4 are used due to the symmetry of the coset. In order to quantify this sensitivity we 
study directly the cross section for these processes in the next section. 



4 Sensitivity of 2 — > 3, 4 cross section to anomalous cou- 
plings 

In this section we analyse the cross section for the 2 — > 3, 4 processes at the parton level for 
couplings using the Lagrangian Eq. ^ implemented in CalcHEP. We studied several different 
channels but will report only on the most representative ones. 

In order to show examples of the enhancements that result from the anomalous Higgs 
couplings, we compute cross sections with a Higgs mass = 125 GeV. We implement a 
cut in the invariant mass of the final state pions m + _ > 200 GeV, such that the Higgs is 
never on-mass-shell when coupled to two pions, as is the case in the SM (the Higgs does not 
decay to a on-shell pair of gauge bosons). We show in Figjljthe ratio of the cross section as 
a function of a (keeping the other parameters fixed in one case and, in the case of MCHM4, 
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changing them according to Eq.(ll )) to the SM cross section at a fixed center-of-mass energy 

of \/s = 2 TeV. Several channels are shown, such as (00, H — ) — >• H — , H — h, hhh, H 1 — . 

The notation (00, H — ) indicates that both 00 and H — initial states have been taken into 
account in the cross section. Thick lines show results for different values of a but keeping 
the other parameters fixed at their SM values whereas thin lines show results when the other 
parameters are changed according to the MCHM4 where, in addition to the relations given in 
Eq.(ll), one assumes d 3 = ^1 — £. Of course not all the values of the anomalous couplings 
are allowed: this plot is only meant for illustrative purposes. 



(00,+-) >+- 

(00,+-)^+-H (00,+-)^ +-H (MCHM4) 




-1 -0.75 -0.5 -0.25 0.25 0.5 0.75 1 



Figure 1: Ratio of the SMEFF (thin lines) and MCHM4 (thick lines) cross sections to the 
SM one versus a parameter at a fixed energy of ^ = 2 TeV. The different channels are: 

(00, H — ) ->• H 1 — (dashed line), (00, H — ) ->• + - h (dot-dashed line), (00, H — ) -> hhhh 

(solid line), and (00, H — ) — > H — (dotted line) for comparison. The notation (00, H — ) 
indicates that both 00 and H — initial states were taken into account. 



Large enhancements of the order of 10 6 with respect to the SM value are easily obtained 
even for small deviations (as small as 10%) of the couplings from their SM values. The 
behaviour of the curves are easily understandable: the 4, 3 and 2 dips in the cross section 
versus a for 2 — > 4, 2 — > 3 and 2 — > 2 are due to the 4 , 3 rd and 2 nd order polynomials 
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in a in the amplitudes. One can see that the enhancements in 2 — y 2 processes are modest 
compared to processes with higher multiplicities, at least at 2 TeV. For the MCHM4 case 
there is a suppression in the 2 — y 3 process, as expected from the parity symmetry of the 
coset. Since the MCHM4 always predict smaller deviations, in what follows we will consider 
the more optimistic case where the parameter a can be the only one different from the SM 
values. 

Next we study the growth with center-of-mass energy of the cross section for different 
multiplicities for a few values of the anomalous coupling, namely a = 0.9, 0.95 and 1 (SM), 
keeping the other couplings at their SM values. We believe that these values of the anomalous 
couplings can be representative of the behaviour of the cross sections, that is, we expect 
the same order-of-magnitude enhancements if the other couplings are also anomalous (but 
without obeying the MCHM4 relations). In Fig. [2] we present a comparison of the cross 
section as a function of energy among representative processes with 2, 3 and 4 particles in 
the final state, for different values of the anomalous coupling. We also show the unitarity 
limit Eq. ([!]). A few comments are in order. The SM cross section quickly stabilises at a small 
value, which depends on the specific process (of the order of 10~ 3 pb, 10~ 2 pb and 10 _1 pb 
for 00 hhh, (00, H — ) — >■ H — h and (00, H — ) —y H — , respectively). It is not surprising 
that in the non-SM case the cross sections grow very fast with energy, reaching up to order 
of 100 pb and violating unitarity at center-of mass energies of the order of a few TeV. It is 
also anticipated that larger multiplicity processes, in the absence of a complete cancellation 
mechanism, grow faster due to phase space. However, what is somewhat unexpected is the 
energy scale at which multiparticle cross sections become comparable to 2 — y 2 processes. 
In the examples shown the 2 — y 3 start to overcome 2 — y 2 at energies of C(lTeV). This 
may be signalling the onset of nonperturbative behaviour well before the unitarity bound is 
reached. It is not clear whether new physics must come in at these scales, as for instance 
the appearance of new resonances. In this work we assume that this is not the case. We also 
checked that the 2 — y 4 process grows very rapidly for a ^ 1, but since it starts out very 
suppressed it surpasses the 2 — y 2 only at very high energies, of the order of 0(5TeV). 

5 Cross sections in the SM with anomalous Higgs cou- 
plings 

So far we have only analysed the scattering of the longitudinally polarised gauge bosons. 
Since it is difficult to separate out the contributions from these polarisations in an experi- 
mental setting, it is important to understand how the large enhancements found will affect 
the corresponding unpolarised cross section. In order to do so, we use the Lagrangian Eq.(|9]) 
promoting the partial derivatives to full covariant derivatives and adopt the unitary gauge 
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(00,+-) -> +- (a=0.9) (00,+-) -¥ +-H (a=0.9) 

(00,+-) ^ +- (a=0.95) (00,+-) -> +-H (a=0.95) 

(00,+-)^ +-(a=1.0) (00,+-)^+-H(a=1.0) 




500 750 1000 1250 1500 1750 2000 2250 2500 2750 



Vs (GeV) 



(00,+-) -> +-+- (a=0.9) (00,+-) -> HHH (a=0.95) 

(00,+-) ^ HHH (a=0.9) (00,+-) -> HHH (a=1 .0) 




500 750 1000 1250 1500 1750 2000 2250 2500 2750 

Vs (GeV) 

Figure 2: Comparison among cross sections as a function of the center-of-mass energy for 
processes with 2,3 and 4 particles in the final state. In the top plot, the solid lines are for 
(00, H — ) — > H — for a = 0.9 (thick), a = 0.95 (medium thick) and a = 1 (thin). Dashed 
lines are for (00, H — ) — > H — h, with the same pattern for the thickness of the lines. In the 
bottom plot, the same pattern of lines show the results (00, H — ) — > hhh and the process 

00 —M 1 — is shown as a dashed line for a = 0.9. In these plots only a deviates from the 

SM value. The unitarity bound is shown as a shaded area in the top right corner. 
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(U = 1). As an illustration, we fix the partonic center-of-mass energy at 2 TeV and com- 
pare the cross sections for the longitudinally polarised gauge bosons with the full SM in 
Table [1} We still keep the notation 0, +, — to indicate longitudinally polarised gauge bosons 
and Z, W ± to denote the unpolarised gauge bosons. 



channel 


a — o — 1 lom.) 


a = U.y; 0=1 


00 +- 


0.13 


295 


ZZ ->■ w + w~ 


610 


655 


00 -> + - h 


2.0 x 10~ 3 


350 


ZZ ->■ W+W~h 


10.9 


46.2 


00 -> /i/i 


0.18 


158 


ZZ -»■ Wi 


7.61 


15.7 


00 -> /i/i/i 


4.9 x 10" 4 


112 


ZZ ->■ h/i/i 


4.65 x IO" 2 


13.6 



Table 1: Comparison of 2 — y 2 and 2 — > 3 cross sections (in picobarns) at y/s = 2 TeV. 



Notice that the processes with longitudinal polarisations are subdominant in the SM. 
However, as we discussed above, they are greatly enhanced with small deviations of the 
couplings and actually dominate the cross sections. The results are consistent with the fact 
that a a \i ~ oll/9 when the contribution from longitudinal polarisations is dominant. The 
enhancements are larger when the final state multiplicity is larger, as expected. For instance, 
cr(00 — y H — h) > cr(00 — y -\ — ) for a = 0.9 and 6 = 1. When all polarisations are included, 
the contribution from transverse polarisations can mask the increase in the cross section for 
the longitudinally polarised gauge bosons. This can be seen in the case of ZZ — y W + W~ , 
where the total increase in the cross section is less than 10%. On the other hand, in cases 
where the contributions from the transverse polarisations are not large, as in the case of 
ZZ —y hhh, enhancements of O(10 3 ) can be obtained. Therefore, multiple Higgs production 
offers the best channels to study anomalous couplings. For the unpolarised case, the cross 
section for ZZ —y W + W~ is still at least one order of magnitude larger than the typical 
2—7-3 processes but ZZ —y hh is of the same order as the ZZ — y hhh cross section. In the 
next Section we discuss the impact of these results for the LHC and future colliders. 
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6 Impact of multiparticle production at the LHC and 
future colliders 



In order to estimate the impact of these enhancements found at the parton level arising from 
anomalous Higgs couplings, we have performed a full calculation of pp — > jj + X, where 
j = u,u,d,d,s,s and X = W+W~ ,W + W~h, hhh at the LHC (for yfi = 14 and 33 TeV) 
using Madgraph5 (vl.4.8). We use CTEQ6L1 parton density function for the evaluation of 
the tree-level cross sections with the QCD scale equal to M z . The selection and acceptance 
cuts include the requirement of two jets with P? > 30 GeV with l^-j < 5 separated with 
AR = A0|_. + Arf-- > 0.4. Besides the cross sections evaluated for the cuts above we have 
evaluated another set of the cross sections for an additional cut to select the vector boson 
fusion process by requiring each jet to be quite energetic with Ej > 300 GeV as well as a 
large rapidity gap between the two jets |Ar/jj| > 4 (see e.g. [9] for detailed motivation of this 
choice). In Table [2] we present the results with and without the vector boson fusion cut. 





14 TeV 


33 TeV 


Process 


with (without) VBF cuts 


with (without) VBF cuts 




a=1.0 


a=0.9 


a=1.0 


a=0.9 




b=1.0 


b=1.0 


b=1.0 


b=1.0 


pp -> jjW + W~ 


95.2 


99.3 


512 


540 




(1820) 


(1700) 


(5120) 


(5790) 


PP -> jjW + W~h 


0.011 


0.0088 


0.0765 


0.0626 




(0.206) 


(0.172) 


(0.914) 


(0.758) 


PP J 3 hhh 


1.16 x 10~ 4 


0.0566 


0.00115 


1.85 




(3.01 x 10" 4 ) 


(0.0613) 


(0.00165) 


(1.46) 



Table 2: Cross section (in fb) for pp — > jjW + W , pp — >• jjW + W h and pp — > jjhhh 
processes evaluated with Madgraph5. 



One can notice that when there are gauge bosons in the final state the cross section 
actually decreases for most cases with a = 1 versus a = 0.9 ones. This is because we chose in 
our example a < a SM = 1 and since the transverse polarisations dominate the cross section, 
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reducing the coupling a results in a smaller cross section. However, in the case of triple 
Higgs production, the enhancements are substantial: roughly a factor of 500 for y'i = 14 
TeV (LHC14) and 1600 for = 33 TeV (LHC33), with VBF cuts. We show in Fig. g the 
results for the triple Higgs production cross section for both LHC14 and LHC33 for the large 
range of anomalous coupling a. The enhancements with respect to the SM case a = 1 are 
large and don't change significantly once \Aa/a\ > 0.1. 

jjhhh with VBF cuts 

33 Tev, MCHM4 33Tev 




-5 : 

10 r 

0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 

a 



Figure 3: Cross section for triple Higgs production pp — > jjhhh with VBF cuts as a function 
of the anomalous coupling a for LHC14 (dark lines) and LHC33 (light lines). Dashed lines 
are for other parameters fixed to SM values and solid lines are for parameters given by 
MCHM4 relations. 

Though the enhancement can be as large as 10 5 for a = 1.5, the absolute value of the 
respective cross sections are quite low (about 10 fb for yfs = 14 TeV with VBF cuts) making 
the study of these processes challenging at the LHC. A dedicated analysis (which are outside 
the scope of this paper) are necessary to understand the LHC or LHC33 sensitivity to the 
processes above. However we can already estimate that if such sensitivity is possible it can 
take place only at high luminocities and/or high energies for quite large values of the a 
parameter. A detailed analysis of these processes at future e + e~ colliders is being performed 
in [19]. 
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7 Conclusion 



In this Letter we have studied multiparticle production in models with anomalous Higgs 
couplings, such as the composite Higgs models. The modified couplings result in a partial 
unitarization of the scattering amplitudes. We found that, due to the stronger cancellations 
in the corresponding amplitudes compared to the usual 2 — > 2 processes, very large enhance- 
ments with respect to the SM, as large as O(10 6 ), can arise at the parton level in the cross 
section of longitudinally polarised gauge bosons, even for small deviations of the couplings 
from their SM values. The cross sections grow faster with energy for larger multiplicities, as 
expected from naive phase space considerations. We pointed out that some 2—7-3 processes 
become as important as 2 — > 2 processes for relatively low energies, of the order of a TeV, 
signalling the onset of nonperturbative effects. When accounting for the contributions from 
the transverse polarisations, the enhancements are somewhat diluted but remain important 
in some processes, especially triple Higgs production. 

However, we showed with a realistic calculation that even with these large enhancements 
the search for multiparticle processes will remain a challenge for the foreseeable future. On 
the other hand, multiple gauge and Higgs boson production receives large enhancements 
in the case of anomalous Higgs boson couplings, and its study could be an important part 
of future experimental programs aiming at understanding underlying theory beyond the 
Standard Model. 
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